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SUMMARY 



A comparison is made of free-molecule— flow theory with the results 
of wind-tunnel tests performed to determine the drag and temperature— * 
rise characteristics of a transverse circular cylinder. 1 

The model consisted of a 0.0031— inch-diameter cylinder constructed 
of butt-welded iron and constantan wires. Temperatures were measured at 
the center and both ends of the model. The cylinder was tested under 
conditions where the Knudsen number of the flow (ratio of mean— free- 
molecular path to cylinder radius) varied between 4 and 185 . The molec— . 
ular speed ratio (ratio of stream speed to the most probable molecular 
speed) varied from 0.5 to 2.3. In terms of conventional flow parameters, 
these test conditions corresponded to a Mach number range from 0.55' to 
2.75 and a Reynolds number range from 0.005 to 0.90* 

’ v The measured values of the cylinder center-point temperature con- 
firmed the salient point of the heat— transfer analysis which was the 
prediction that an insulated cylinder would attain a temperature higher 
than the stagnation temperature of the stream. 

. Good agreement was obtained between the theoretical and the experi- 
mental values for the drag coefficient. As predicted, theoretically, 
the drag coefficient was independent of Knudsen (or Reynolds number) 
variation and depended primarily upon the molecular— speed ratio. 


INTRODUCTION 


The various phenomena connected with the high-speed flows of rare- 
fied gases have- recently been the subject of several analytic investi- 
gations (references 1, 2 , 3 , and 4). Interest in this phase of aero- 
dynamics and thermodynamics has been aroused primarily because of the 
necessity of being able to compute forces on and heat transfer to. mis- 
siles which may fly at great altitudes at high speed. In another 
analytic exploration of this field, simplified- concepts of aerodynamic 
and thermodynamic phenomena occurring at high altitudes have been used 
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"by Whipple (reference 5) in an investigation of the properties of N the 
upper atmosphere "by means of measurements of the brightness and deceler- 
ation of meteors. 

Experimental work in the field of high-speed rarefied gas flows has 
been practically nonexistent up to this point. The pioneer work of 
Epstein (reference 6) on the drag of spheres in a rarefied gas was con- 
cerned only with very low— speed motions. Consequently, the purpose of 
the present investigation was "twofold: (l) to obtain exper ime ntal data 

that could be compared with the predictions of an analytic treatment of 
aerodynamic and thermodynamic processes occurring at very high altitudes; 
and (2) as a means to this end, to develop a wind-tunnel and related 
instrumentation that would permit testing of models under simulated con- 
ditions of high— altitude high-speed flight. In connection with the second 
objective — the development of a wind tunnel and instrumentation — it 
should be pointed out that a low— density wind tunnel is being developed 
concurrently at the University of California and the authors wish to 
acknowledge the assistance they have received from the. reports issued by 
this group, as well as the stimulus received from conversations with the 
University of California personnel. 

At the beginning of the experimental work, it was decided that the 
initial model tests in the low-density wind tunnel should be confined to 
the free-molecule— flow regime despite the fact that this regime is pro- 
bably of less technical importance than the slip— flow and intermediate 
regimes. This decision was reached because the free-molecule regime is 
susceptible to analysis; whereas the slip— flow and intermediate regimes 
have not yet been adequately treated (cf . , reference 2 ). In view of 
the many uncertainties connected with the development of entirely new 
types of wind-tunnel instrumentation, it was thought to be essential to 
test in a field where the experimental work could be guided by a theo- 
retical background. 

The selection of the test body was governed both by the requirement 
that free-molecule flow be obtained and by the practical considerations 
arising from tunnel size and instrumentation. The former condition re- 
quired that the significant body dimension be small compared to the 
mean— free-molecular path. 

A cylinder with axis perpendicular to the stream was chosen as the 
body which most nearly fulfilled the requirements of the test. The 
radius of the cylinder could be made small compared with the mean— free 
path. In addition, since it was desired to measure * both aerodynamic 
forces and body temperature, the cylinder had the virtue of being doubly 
symmetrical so that only drag force had to be measured. The problem of 
temperature measurement could be solved by constructing the model from 
dissimilar metals so that the model itself formed a thermocouple. 

The present investigation is a first step toward an attempt to 
gain an understanding of the comparatively new and unexplored' field of 
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high-speed loir-density flows. ,The experiments reported herein were made 
to provide a comparison with some important results of free-molecule— 
flow theory and to facilitate the solution of problems connected with 
the design, operation, and instrumentation of a low-density wind tunnel. 
The tests were conducted in the Ames low— density wind tunnel. 

NOTATION 


surface area, square feet 


c x 

Cy 

c z 

Cp 

pDi 

CD r 

D 

Ei 

E r 

E r 

Et 

f(s) 

f°(s) 

g(s) 

g°(s) 

G 

lo 

*1 

J 

k 


components of total velocity, feet per second 

total drag coefficient, dimensionless 
impinging drag coefficient, dimensionless 
re— emission drag coefficient, dimensionless 
total drag force, pounds 

rate of incident molecular energy, foot-pounds per second 
rate of re— emitted molecular energy, foot-pounds per second 
rotational energy, foot-pounds per molecule 

rate of incident translational molecular energy on front side of 
body, foot-pounds per second 

dimensionless function of s defined by f ( s)=Z 1 (b 2 +3)+Z 2 ^ s2+ i) 
dimensionless function of s defined by f°(s)=Z 1 (s2+2)+Z 2 ^ s 2 +|^ 
dimensionless function of s defined by g(s) = 3(Z 1 +Z 2 ) 
dimensionless function of s defined by g°(s) = 2(Z 1 +Z 2 ) 
normal force, pounds 

modified Bessel function of first kind and zero order 
modified Bessel function of first kind and first order 
mechanical equivalent of heat, 778 foot-pounds per Btu 
Boltzmann constant, 5 .66 x 10 24 foot-pounds per °F per molecule 
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Knudsen number (ratio of meanr-free— molecular path to cylinder 
radius) , dimensionless 

direction cosines, dimensionless 

length, feet 

mass of one molecule, slugs 
Mach number, dimensionless < 

number of molecules striking body, molecules per square foot, 
second 

number of molecules per unit volume of gas 

free— stream static pressure, pounds per square foot 

/ 

cylinder radius, feet 

universal gas constant, 15^4 foot-pounds per pound— mole , 

°F absolute 

molecular- speed ratio, (ratio of stream mass velocity to most 
probable molecular speed) , dimensionless 

free— stream static temperature, °F absolute 

tunnel stagnation temperature, °F absolute 

cylinder temperature , °F absolute 

surface temperature , °F absolute 

test— section vail temperature, °F absolute 

mass velocity, feet per second 

components of mass velocity, feet per second 

re-emission velocity, feet per second 

most probable molecular speed, feet per second 

molecular weight, pounds per pound— mole 
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x 

yi Cartesian .coordinates . 

z 

J 

Z\ dimensionless function of s defined by Zi = ite-® 2 / 2 I 0 (s?/2) 

Z 2 dimensionless function of s defined by 

Z 2 = ns 2 e- s2 / 2 I 0 (s^/2) + I x (s s /2) 

a accommodation coefficient, dimensionless 

3 reciprocal of most probable molecular speed, seconds per foot 

r gamma function 

7 ratio of specific heats, dimensionless 

e emissivity, dimension?, ess 

0 angle of body element with respect to free stream, radians 

p density of free stream, slugs per cubic foot 

0 Stefan— Boltzmann constant, 3*74 x 10~ 10 foot-pounds per sq uar e 

foot, V 

. dimensionless group defined by 

1 + — s sin 0 [l + erf (s sin 0)] e s2 s i Q2 ® 

\lr = 1 + — ' 

o' P p 

1 + vjt s sin 0 [l + erf (s sin 0)] e s G ^ n ® 
dimensionless group defined by 

1 — s sin 0 [ 1 — erf (s sin 0)3 e s2 sin 2 0 

V = i + -4 — : 

1 — vjt s sin 0 [l ~ erf (s sin 0)] e s2 sin 2 0 
X dimensionless group defined by 

X s e™ 62 sin 2 0 + */~Jt s sin 0 [l + erf (s sin 0) 3 

X dimensionless group defined by 

X = e“^ 2 sin 2 0 s gin* q [i — e rf (s sin 0)] 

erf (a) error function ( — P e~ x2 dx J, dimensionless 

\ a,/ i ^ ✓ 
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ANALYSIS 


The general problem of the calculation of aerodynamic forces and 
heat transfer to a body in a free-molecule— flow field is amenable to 
solution by application of the simpler concepts of the kinetic theory 
of gases. The simplification arises from the assumption that the motions 
of molecules that strike the body are in ho way affected by collisions 
with molecules that have struck the body and have returned to the main 
stream of gas. This basic postulate allows the molecular motions of the 
oncoming stream to be treated as having the classical Maxwellian velocity 
distribution superimposed on a uniform mass velocity. Hence, the severe 
complication of dealing with a nonuniform gas is avoided. Of course, 
this assumption does not imply that returning molecules do not collide 
with oncoming, molecules since collisions do occur, on the average, at a 
distance from the body equal to one mean free path. However, if this 
mean free-path length is large compared with the body dimensions, it can 
be seen from geometrical considerations that the probability is small 
that the impinging molecules that do collide with returning molecules 
will strike the body without first having had their original Maxwellian 
velocity restored by collisions with other impinging molecules. It would 
seem probable, also, that a high value of mass velocity compared with the 
random thermal velocity of the gas would favor free-molecule flow because 
the re— emitted molecules would tend to be swept along with the stream 
away from 'the zone of influence of the body. 

In general, the mechanism of energy and momentum transfer for 
free-molecule flow are inseparable. In effect, this means that the drag 
on a body is affected by the surface temperature of the body. This state 
of affairs arises because of the fact that a portion of the momentum, 
imparted to the body is due to the emission ctf molecules from the sur- 
face and the velocity of emission depends upon the temperature of the 
body. This phenomenon will be discussed in detail later. 

Because of the assumption that the motions of the oncoming and 
re-emitted molecules can be treated independently, the computation of 
the drag on and heat transfer to a cylinder can be broken down into 
separate parts — the part due to impact by impinging molecules and a 
part due to re-emission of the molecules from the surface. By the use 
of this basic principle, it is possible to calculate the drag and heat- 
transfer characteristics of a transverse cylinder in a free-molecule— flow 
field. 


Heat Transfer to a Transverse Cylinder 
in a Free-Molecule-Flow Field 


The calculation of the heat— transfer or energy— exchange process 
between the cylinder and the gas is made by using the basic method 
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described in reference 4. The general equations of- reference 4 are 
expressed in Cartesian coordinate form and can be summarized as follows: 

The translational molecular energy incident upon the front side of 
an elemental plane area dA inclined at an angle 6 with respect to the 
stream mass velocity U is, for a monatomic gas, 

dEt = n dA (l) 


where n, the number of molecules striking unit surface area per second, 
is given by 


n 


Nv m x 
2 V ji 


( 2 ) 


For the rear side of the elemental area, similar expressions apply. 


dE t ' 


mir 

2 


+ y'kT ) dA 


(3) 


The total molecular energy incident upon the front and rear sides 
of dA is, for a monatomic gas, given hy equations (l) and (3), respec- 
tively, since only translational energies are involved. However, for a 
diatomic gas, the expression for the incident energy contains an additional 
term due to the rotational, component of the molecular internal energy, and 
this component must he added to the translational component. The rotat- 
ional energy per molecule of a gas composed of rigid dumbbell molecules is 


E R = kT 


(5) 


Then, the total incident molecular energy on the front side of dA 
is, for a diatomic gas, given hy the following expression: 


dEi = dEt + oEfidA. 



(6) 


and, for the rear side of dA, 


dEi’ = dE t * + n*E R dA = n' 


+ 0|r’+l) kT 


dA 


(7) 


The method of calculation of the energy transported from the elemental 
area hy re-emitted molecules is shown in appendix A. The result derived 
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"therein can he applied to the calculation of the energy transport hy a 
monatomic gas from the front side of the area hy combining equations (Al) 
and (A12) to yield 

dE r = dEt(l-a) + 2ankT c dA (8) 

and, for the energy transport hy a diatomic gas, as 

dE r = (dEt + nkT dA) (l-a) + 3ankT c dA (9) 

Parallel equations can he written for the rear side of the area, dA. 

The net radiant energy exchange is- (assuming gray-body radiation) 


ea (T c 4 - T t 4 


dA 


( 10 ) 


where T^ is the effective temperature of the surroundings to which the 
area dA is radiating. This expression implies that the body of which 
dA is a part is small in comparison with the surrounding volume. 

The problem now becomes one of applying the foregoing equations, 
which were derived for a plane area placed at an arbitrary angle with 
respect to the stream, to a specific body — in this case a cylinder with 
axis. perpendicular to the stream velocity. For a cylinder of length L, 
the elemental area may be written as 


dA = rLd© 


(ID 


Then, if- internal energy input to the cylinder per unit area is denoted 
as Q (in the usual experimental arrangement this would be electrical 
energy) , an equation can be written which expresses an energy balance on 
'the cylinder. For a monatomic gas, the energy-balance equation is 


2arL 


f / n (^2~ + * kT ) 


de + 


iTT/2 


+ t'kT^ 


d0 


4arLkT c ^ nd0 + 


r 


n’cLG) + 2nrL 


Q - ea ^T c 4 - T-t 4 \ 


= 0 


( 12 ) 
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For a diatomic gas, a similar equation applies 




The reduction of the integrals appearing in equations (12) and (13) is a 
straightforward tut lengthy procedure . The methods used to evaluate these 
integrals are given in appendix B. The final results are written here, 
for a monatomic gas,. as 


2 — (Z 1+ Z a ) - [z^s^) + z / b s + |)1 + [ ec(T c 4 -Tt 4 ) " Ql = 0 

L \ 4 /J l_ J ( 14 ) 

and, for a diatomic gas, as 



In order to facilitate computation, values of the dimensionless groups, 

f(s)= Z!(s 2 +3)+Z 2 ^b 2 + 1 ^, g(s)= 3 (Zi+Z 2 ), f°(s)= Z x ( s 2 +2)+Z 2 ^s 2 + |^, 

and g°(s) = 2 (Zj+Z 2 ), appearing in equations (14) and (15), are tabulated 
in the table of appendix B. 


Equations (l4) and (15) then completely describe energy transfer 
occurring between a transverse cylinder and its surroundings in terms of 
the cylinder temperature', static gas temperature and pressure , molecular- 
speed ratio, accommodation coefficient, the effective temperature of the 
surroundings to which the cylinder is emitting radiation, the emissivity 
of the cylinder surface, and the internal energy input to the cylinder. 

The equations are subject to the following restrictions: 

1. The surface temperature of the cylinder is assumed to be constant, 
circumferentially. 

2. The cylinder receives and emits radiant energy as 1 a gray body; 
that is, the emissivity is independent of wave length. Also, the cylinder 
is assumed to be small compared to the surrounding volume . 

3. The internal energy of the gas stream in which the cylinder is 
located has the value (3/2)kT for a monatomic gas and (5/2)kT for a 
diatomic gas . 
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Drag on a Transverse Cylinder in a Free-Molecule— Flow Field 


Hie general problem of aerodynamic forces acting on bodies in free— 
molecule flow has been treated by several authors (references 1 , 2 , 
and 3). Hie basic method of calculation was identical in each case. 

Hie method consists of the assumption of a gas with a Maxwellian velocity 
distribution superimposed upon' an arbitrary mass velocity. Hie number 
of molecules having specified velocities that impinge upon an elemental 
area inclined at* an arbitrary angle to the stream mass velocity is com- 
puted and , from this calculation, the force due to the transmitted 
momentum is obtained* To this momentum due to molecular impingement, 
however, must be added the momentum imparted to the surface by molecules 
emitted from the surface. The magnitude and. direction of the force 
imparted by emitted molecules depends, among other things, on the type 
of emission that is assumed to occur. A condensed discussion of the 
several possible types of surface emission can be found in reference 4. 

An extensive treatment of the subject is given in reference 7. 

In order to determine the total aerodynamic force on a body of 
arbitrary shape, the differential momentum imparted to an element of the 
body surface is computed, as outlined in the preceding paragraph, and the 
total momentum obtained by integration over the body surface. Also, since 
the impinging and re— emitted momenta can be added linearly, the resulting 
total momentum can be calculated for any desired type of re-emission. 

■ In this paper, completely diffuse scattering of the incident molec- 

ular stream from an insulated surface is assumed. By this assumption it 
is meant that (a) the molecules that are scattered from the surface are 
emitted with a Maxwellian distribution of speed, (b) the direction of 
molecular emission from the surface is controlled by the Knudsen cosine 
law (cf., references 7 and 8), (c) the body temperature is assumed to be 
uniform over the surface , and further (d) that conservation of number of 
molecules is maintained for each element of the body surface. These 
assumptions seem to be in accord with physical fact, as discussed in 
reference 6. 

The specific problem of the aerodynamic drag on a transverse cyl- 
inder has been treated in references 1 and 3 * The reason for repetition 
of the calculation here is that the results have been obtained in a 
compact closed form which is amenable to easy calculation and, further, 
that a solution for the re— emission drag is obtained which involves 
measurable physical quantities. Appendixes A and C contain the details - 
of the calculations and only the final results are presented here. Hie 
results of the> conputations are as follows: The total drag acting on 

an insulated transverse circular cylinder of radius r and length L 
is given by 

D = (2rL)(C Di + Cp^) 


(16) 
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where and Cp r are the drag coefficients associated with impinging 

and re— emitted momenta. The value of is shown in appendix C to 

have the value 

CD* = — -jl 0 (s 2 /2) + (^Y~) [lo(s 2 /2) + I x (s 2 /2) | (17) 

and the value of Cj> r for an insulated cylinder is shown in appendix A 
to he, for a monatomic gas. 



and, for a diatomic gas, 

(19) 


APPARATUS AND TEST METHODS 


The Wind Tunnel 


The Ames low— density .wind tunnel , in which the present tests were 
conducted, is an open— jet nonreturn— type tunnel, A sketch showing /the 
major components of the wind tunnel is shown in figure 1. The test gas, 
from a bottled source, is throttled into a settling chamber and then 
passes through an axially symmetric nozzle which discharges a free jet 
into an 18— inch cubical test chamber. The test chamber is connected 
through a surge chamber to four boos ter—' type oil-diffusion pumps which 
continuously evacuate the system. The diffusion pumps are connected in 
parallel and axe backed by a mechanical vacuum pump which discharges to 
the atmosphere. Valves are provided at the inlet to each diffusion pump 
in order to allow the pumping system to be isolated when the test chamber 
is opened to the atmosphere. 

A traversing mechanism is installed in the test chamber which permits 
three degrees of motion relative to the gas stream. The mechanism is used 
to position models and instrumentation in the desired location within the 
gas stream. 

The nozzles used in the tests were axially symmetric with an outlet 
diameter of 2 inches and were machined from clear plastic. The nozzle 
contours were determined by the method outlined in reference 9. No 
boundarjr-layer correction was applied to the nozzle contours. 


jjS/2 /z 1 (s 2 +3) + Z ; 


bs N 


>(**+ i) 


3(Zi+z 2 ) 
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After the system vas evacuated, the flov rate of the test gas into 
the settling chamber was adjusted until the stagnation pressure, as 
measured by a McLeod gage, reached the desired value. 

The static pressure in the test chamber vas read on another McLeod 
gage and, if the resulting pressure ratio was that required- to establish 
flow at the particular Mach number desired for the run, the tunnel was 
considered ready for tests. In order to match tunnel conditions of 
static pressure and pressure ratio with those that existed at an earlier 
time, it was sometimes necessary to alter the amount of throttling at the 
pump inlets due to a slight -variation of pumping speed from day to day. , 

, The pressure level in the tunnel was varied by increasing or 
decreasing the flow of the test gas into the system. The static and 
stagnation pressures were again checked to determine if sufficient pres- 
sure ratio was available to establish flow at the desired Mach n umb er. 
This procedure was somewhat complicated by the fact that at the lower 
pressure levels the boundary layer in the nozzle would tend to thicken^ 
thus causing a decrease in the Mach number compared with that obtained 
at the higher pressure level's. 

The wind tunnel could be operated through a range of test— section 1 
pressures of from 20 to 135 microns of mercury for a given nozzle and 
over a range of Mach numbers from 0.5 to 2.75 by using several nozzles. 


Instrumentation and the Test Model 


The static pressures which existed in the wind tunnel were, in 
general, too low to be measured with a conventional manometer or Bourdon- 
type pressure gage; therefore, two types of gages widely used in high- 
vacuum work were utilized, namely, the McLeod and Pirani gages. For a 
complete description of these gages the reader is referred to any 
standard text on high-vacuum technique (e.g., reference 10); however, 
some of the special features of the gages used will be described. 

The McLeod gage used as a primary standard in these tests had an 
initial volume of l46 milliliters connected to a closed tube of 1.95- 
millimeter diameter. This arrangement gave a scale length of 
centimeters for a pressure range Of 0 to 3000 microns. (For the sake of 
brevity, the term micron is used throughout as a unit of pressure rather 
than the term microns of mercury. ) The diameter of the closed tube was 
measured by the mercury pellet method to within ±0.01 millimeter and the 
difference in height of the mercury columns could be read to the nearest 
0.1 millimeter by use of a micrometer scale. At a pressure of 100 
microns, this gage could be read to within 1 percent with ease. 
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The gage was connected to the test chamber with an 18-mill i meter— 
diameter glass tube to insure a short time lag between equalization of 
the system and the gage pressures. A dry— ice and acetone cold trap was 
provided between the gage and the test chamber to prevent condensable 
vapors from reaching the gage. 

In general, it is difficult to obtain a completely independent 
method of checking the accuracy of pressure measurements made with a 
McLeod gage for pressures below about 1 millimeter of mercury. Above 
1 millimeter of mercury a U— tube manometer, f ill ed with a very .low— vapor- 
pressure oil and having a high— vacuum system (of the order of 10 — 5 mm 
of Hg) connected to one leg, may be used. The McLeod gage used in the 
tests was checked in this manner and the agreement between the oil 
manometer and the McLeod gage was within ±3 percent.' 

An indication of the accuracy of the McLeod gage used in the tests 
in the pressure range below 200 microns was obtained by comparing the 
pressures indicated by the test McLeod gage with those obtained with a 
laboratory— type McLeod gage having a range of 0 to 200 microns, which 
has been calibrated against 'a laboratory standard McLeod gage at the 
manufacturer’s plant. The agreement between the test gage and the 
laboratory— type gage was within ±3 percent at a pressure of 100 microns. 

The Pirani gage vised in the tests, consisted of an open glass tube, 
which contained a heated filament. The gas pressure in the tube deter- 
mined the temperature and thus the resistance of the filament. A 
temperature— compensating tube of identical construction was sealed off 
at a very low pressure. The measuring and compensating tube filaments 
formed two legs of a Wheatstone bridge; and the unbalance of the 
Wheatstone bridge circuit, as indicated by a self-balancing potentiom- 
eter, wa,s a measure of the pressure in the open leg of the gage. The 
Pirani gage was calibrated against the 0— to 3000-micron— range McLeod 
gage as a standard. 

A drift of the calibration of the Pirani gage with time was noted 
and for this reason the Pirani gage was calibrated before and after every 
test run. A typical calibration curve is shown in figure 2. 

Examination ' of figure 2 reveals that the sensitivity of the Pirani 
gage used was nearly constant from a pressure of 10 microns to a pressure 
of 80 microns. Above about 200 microns this Pirani gage becomes too 
insensitive to be usable as a pressure-measuring instrument. 

Whenever a pressure-measuring device is connected to a vacuum system 
by a tube or small orifice, two troublesome effects are present — time 
lag and gassing. The time lag of the pressure gage is due to the fact 
than, as the pressure in the main part of the system is changed, a volume 
of gas must flow out of or into the pressure gage and its connecting 
tubing. At the low pressures encountered during the tests, this time 
lag may be of the order of hours if the gage volume and connecting tubing 
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are not properly proportioned. Gassing, that is, gas being given 15 ) or 
adsorbed by the surface of the pressure gage and its co nn ecti ng tubing, 
has the effect of causing the pressure gage to indicate a different 
pressure than that vbich exists in the main part of the system. Out- 
gassing effects can be reduced materially by proper cleaning of the 
surface and by beating the surface under a vacuum to drive off water 
vapor and other materials which cling to ordinary surfaces. 

Schaaf, Mann, and Cyr (references 11 and 12) have investigated this 
problem in considerable detail, and show that, if the assumption is made 
that gassing is a function of surface area only, there exists an optimum 
geometry of gage volume and connecting tubing which will result in min i— 
mum time lags. References 'll and 12 were used as guides in the design 
of the impact and static— pressure— measuring devices which were used in 
the tunnel calibration tests, and, in addition, the time lag of these 
pre s sure-me asur ing systems was measured. 

The impact pressures encountered during the wind-tunnel-calibration 
tests were too high to be measured with a Pirani gage and, because of the 
difficulty of obtaining a suitable flexible tube which would have a small 
time lag, the McLeod gage could not be used. These impact pressures, 
therefore, were measured with a glass U-tube manometer filled with a 
low— vapor— pressure oil. The manometer was located inside the test cham- 
ber with one leg of the manometer open and the other end fastened with a 
glass connection to a l/ 8 -inch-diameter spherical-head impact tube. 

Figure 3 shows a diagrammatic sketch of this arrangement. The difference 
in oil level between the two legs, as read by a coordinate-type cathetom- 
eter located outside the tunnel, was the difference between the test- 
chamber static pressure and the impact pressure. The tesi>-chamber static 
pressure was measured by the 0- to 3000-inicron McLeod gage described 
earlier. 

The time lag of the system shown in figure 3 was measured by setting 
the pressure in the test chamber at approximately 3 millimeters of 
mercury and then evacuating the chamber to approximately 10 microns as 
quickly as possible. The time required for the oil level in the two • 
legs to equalize was approximately 15 seconds. This indicates a very 
short time lag for the practical case as the rate of pressure change 
encountered in the wind-tunnel . calibration was small conpared to that 
measured during the time-lag test. The procedure used in the wind— 
tunnel— calibration tests was to continue to take readings of the oil 
level until steady state was established. 

The static pressure of the gas stream from the nozzle was measured 
with a 0.10-inch-diameter ogival probe connected directly to a Pirani 
gage as shown in figure 4. Three 0.0225-inch-diameter holes spaced 120° 
apart were located at a point on the ogive where the pressure is equal 
to free— stream pressure as calculated by the method given in reference 
13. The position along the probe where the surface pressure is equal 
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to free— stream pressure varies somewhat with. Mach, number; however, this 
variation is not. large and static pressures measured by this probe 
agreed with pressures measured with a 4°cone. 

No boundary— layer correction was applied to the ogive, or cone con- 
tour ; consequently, it is possible that some error in the measurement of 
static pressure existed. Figure 5 shows the main dimensions of the 
static— pressure probe used. 

In order that the Pirani gage shown in figure 4 could be calibrated 
against the McLeod gage, which was connected to the test chamber, a 4— 
millimeter-bore stopcock was glass— welded to the Pirani gage to allow 
calibration with a relatively large opening and attendant short time lag. 
As the stream static pressures were measured with the stopcock closed, 
the following tests were conducted to measure the order of magnitude of 
the. time lag of the system. The test chamber was pressurized to 165 
microns and suddenly evacuated to approximately 5 microns. The rate of 
pressure change as. measured by the Pirani gage with the stopcock open 
and with the stopcock closed is shown in figure. 6. The procedure was 
repeated with a rising pressure and the results are also shown in 
figure 6. 

It can be seen from figure 6 that the two curves, representing the 
runs with the stopcock open and closed, are essentially congruent after 
70 seconds for decreasing pressure and after 50 seconds for the increas- 
ing pressure runs. With the stopcock open, the time lag of the Pirani 
gage should be very small (references 11 and 12) and the curve of pressure 
as a function of time includes the time required for the pumping equipment 
to evacuate the test chamber , the time lag in the self-balancing potenti- 
ometer, and the thermal lag in the Pirani filaments. The curve of 
pressure as a function of time for the, stopcock closed includes* all these 
factors plus the additional time lag inqposed by the small holes ' in the 
static probe. The pressure variation during the time-lag tests was much 
more severe than any encountered during the static survey tests; therefore, 
it can be concluded that, by waiting for at least 2 minutes before record- 
ing data, all errors due to time lag would be eliminated. 

The test model .- The model tested was a 0 . 0031-inch--diameter right- 
circular cylinder fabricated from iron and constantan wires and had an 
exposed length of 15/32 inch. The iron and constantan wires were electri- 
cally butt-^welded together to form a thermocouple junction. Figure 7 
shows a microphotograph of the model. It can be seen that the welding 
operation did not appreciably alter the shape of the model. The butt- 
welded junction was located in the center of the model and a 0 . 00 31 — inch 
iron wire was silver soldered to the constantan side of the model at a 
distance of 5/l 6 inch from the center junction. A 0.0031-inch-diameter 
constantan wire was similarly fastened to the iron side of the model. 

These two auxiliary thermocouple junctions were used to measure the 
tenperatu re at the ends of the model. 
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The entire model was covered with fine— grain soot obtained from an 
acetylene flame in order that the emissivity would be constant over the 
total model length and so that an established value of emissivity of 
0.95 could be used to evaluate the heat loss by radiation. 

Thermocouple voltages from the center and end thermocouples were 
read on a null— type laboratory potentiometer having a least count of 
0.01 millivolt. 

The drag balance .— A specially constructed microbalance shown in 
figures 8(a) and 8(b) was used to measure the drag of the model tested. 

The balance consisted of a beam mounted upon an agate knife-edge which 
rested in a steel V-block, as shown diagrammatically in figure 9. The 
model was mounted between two arms extending down from either side of 
the beam. Electrical leads were transferred from the moving beam to the 
balance frame through mercury pools. 

A drag force on the model caused an imbalance of the beam which was 
counteracted by a magnetic— restoring force produced by an armature located 
in a magnetic field produced by an alternating-current coil. An optical 
lever approximately 2 feet long was used to determine the zero position 
of the balance. 

The magnitude of the drag force in terms of magnetizing current was 
determined by hanging known weights on a lever arm attached to the balance 
beam and measuring the current necessary to restore the balance to zero 
position. Figure 10 shows a typical calibration of the drag b alanc e. 

*■ The model was mounted between the extension arms of the drag b alan ce 
and was shielded over all but the center 15/32 inch, as shown in figures 
8(a) and. 8(b). A check was made on the effectiveness of -the shielding on 
the supporting members. A small shield was mounted to the center main 
shield (extension arms) covering the center portion of the model (normally 
exposed) . The completely shielded model was then placed in the gas stream. 
No drag forces were detected. 

It was found that the balance exhibited an erratic calibration shift 
from day to day, which necessitated calibration of the balance before and 
after each drag determination. 

' . Calibration of the Wind Tunnel 

f 

The purpose of the calibration of the wind tunnel was to determine 
the Mach number and static— pressure distributions of the stream in the 
plane at which the model was located - l/k inch downstream of the exit 
plane of the various nozzles used. 

The method used was to survey the stream with the impact- and 
static— pressure probes described earlier. These impact— and static- 
pressure probes, together with their pressure-measuring gages, were 
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mounted 2 inches apart upon a common support from the traversing mecha- 
nism and readings were taken at l/l6-inch intervals from the center to 
the edge of the stream. This method of mounting allowed both the impact 
and static pressures to be measured during the course of one run, thus 
obviating the necessity of matching tunnel conditions oh separate runs. 

Table I lists the nozzles tested, their design Mach number, and the 
tunnel conditions for which the surveys were made. 

The static pressure varied across the stream issuing from each of 
the nozzles tested. For most of the nozzles, the pressure was somewhat 
higher near the center portion of the stream than was the pressure of 
the quiescent gas in the test chamber. Certain of these nozzles, 
however, did not show this general trend. In these nozzles the pressure 
first decreased, then increased, as the probe was moved from the noz zl e 
wall to the center of the • stream. Figures 11 and 12 show this effect. 

In general, the static-pressure variation across the stream was not 
large. The average change over the center l/2 inch of the stream 
amounted to 2 percent, while the maximum deviation was only 5 percent. 

The impact-pressure surveys of all the nozzles had one common 
characteristic; that is, the impact pressure increased from a mini mum 
value at the edge of the nozzle wall to a maximum value hear the center 
of the stream. Hie area of the stream which had a constant, or nearly 
constant, impact pressure varied somewhat from nozzle to nozzle, and 
also to some extent with pressure level. The average variation of impact 
pressure over the center l/2 inch of the stream for all nozzles tested 
was 7 percent with a maximum variation in some cases of 14 percent. 

Figures 11 and 12 show the variation of impact pressure across the stream 
for one nozzle using two gases, nitrogen and helium. 

The Mach number of a gas stream can be computed from measured values 
of the stagnation and static pressure of the stream, provided that the 
flow is isentropic between the points of pressure measurement. However, 
if the static and inpact pressures are measured at the same point in the 
stream, no assumption regarding the nature of the flow need be made in 
order to compute the local stream Mach number. For an inviscid flow, 
the Mach number can be calculated from Rayleigh's equation, the use of 
which involves the assumption of a norma l shock wave ahead of the impact- 
pressure tube and isentropic compression of the gas in the subsonic flow 
behind the shock wave to the stagnation point on the impact-pressure tube. 
In a low— density flow, however, viscous effects are present in the flow 
about a stagnation region, thus invalidating the assumption of isentropic 
compression behind the normal shock wave. Chambre (reference l4) computed 
the viscous correction to be applied to sphericai-head impact tubes in a 
low-density flow, and figure 13 shows the Mach number distribution of 
nozzle 1 as computed by the Rayleigh equation (inviscid) and as 
computed by Chambre' s equation. The Mach number calculated from the 
stagnation and local static pressures were, in all cases, higher than 
those computed using local values of impact and static pressure, thus 
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indicating that the flov through the nozzle was nonisentropic due to 
viscous effects. Chambre’s equation was used to compute Mach number from 
measured static and impact pressures for the tests described herein. 

The variation of impact and static pressure across the stream 
resulted in a stream in which the Mach number was not constant. Figures 
13 and lb show this effect and are generally representative of all noz- 
zles tested. If only the center l/2 inch of the. nozzles is considered, 
the average variation in Mach number amounted to 5 percent; however, at 
certain pressure levels in some of the nozzles the variation was as much 
as lb percent. 

In order to determine if shock waves were present in the center 
l/2 inch of the gas stream, the nitrogen-glow— discharge technique of flow 
visualization was used to examine the flow qualitatively. This technique 
was originally described in reference 15 . The method used to generate 
the glow discharge was to pass the test gas through a screen which was 
connected to a pulsed direct— current power supply. The screen was located 
at the entrance plane to the nozzle and was pulsed to a negative peak 
potential of 1000 volts.. The pulsation frequency was 1000 cycles per 
second. ' 

Figure 15 shows the gas stream from nozzle 1 at a static pressure 
of lib microns. A 5° half-angle wedge was placed in the stream to 
generate a relatively weak shock wave which could be used as a yard- 
stick to evaluate the strength of the nozzle shocks, if present. No 
shock waves were observed in the gas stream at this pressure level if 
the wedge was removed. As the static-pressure level was lowered, the 
shock wave from the wedge became too indistinct to be photographed. 
Observation of the nitrogen glow did not indicate the presence of strong 
shock waves in the nozzles. 

Drag and Heat— Transfer Tests 


The test procedure used in the drag and heat— transfer experiments 
was to calibrate the drag balance with the tunnel shut down and then to 
set tunnel conditions corresponding to those obtained during tunnel- 
calibration tests by matching the stagnation and test chamber static 
pressures. As extremely close control over tunnel conditions was diffi- 
cult to attain, the pressure level' during certain test runs varied 
approximately 5 percent from those obtained during the tunnel— calibration 
tests . Examination of the results of the tunnel— calibration tests showed 
that Mach number variation is small for small changes in the pressure 
level; therefore, the Mach number obtained during nozzle— calibration tests 
was applied to the heat-transfer and drag tests without correction. The 
Mach number used in plotting the data ( which was convertednumerically 
to molecular-speed ratio by means of the relation s * M*/// 2 ) vas tbe 
average value obtained from a graphical integration of the Mach number 
distribution across the center 15/32 inch of the stream. 



NACA TN 2244 


19 


After stable operating conditions had been attained, the drag force, 
model temperatures, and other pertinent data were recorded. Table II 
lists the recorded data. It will be noted that, in the column listing 
the recorded drag force, two values occur on some of the runs. The dif- 
ferent values were due to separate drag measurements during the course 
of one run or resulted from small shifts in the calibration of the drag 
balance as determined before and after each run. In the figures which 
are later presented, showing drag coefficients, the points correspond to 
those calculated using the average of the two values listed in table II. 

To determine if the model surface conditions were changing with time, 
the model drag and temperature were measured at 1 /2-hour intervals over 
a period of 6 hours for one particular test condition. Figures l6 and 17 
show the results of these tests and it can be seen that the drag coeffi- 
cient and temperature ratio remained essentially constant, thus indicating 
that the accommodation coefficient and tunnel conditions were constant 
over this period of time. 

. The reproducibility of the data was checked by repeating the above- 
described run on separate days. The reproducibility was found to be 
within the accuracy of the drag and temperature measurements described 
above . 


RESULTS AND DISCUSSION 


A Discussion of the Assumptions and Results 
of the Theoretical Analysis 


An examination of the heat— transfer equations (l4) and (15) for a 
monatomic and a- diatomic gas, respectively, reveals several interesting 
features. If an adiabatic case is considered; that is, if both internal- 
energy input to the cylinder and radiant— energy exchange are considered 
absent, it can be seen that the ratio of cylinder temperature to free— 
stream static temperature is a unique function of the molecular- speed 
ratio. The values of the temperature ratio- computed from equations (l4) 
and (15) are shown in figure l8. For comparison, curves of the tempera- 
ture ratio for continuum flow are also shown. These latter curves were 
computed for the case where the total stream energy is recovered — in 
other words, the temperature recovery factor was considered to be unity. 
Several remarkable conclusions can be drawn from an examination of this 
figure. It can be seen that an insulated cylinder in a free-molecule- 
flow field will attain a higher temperature at a given value of molecular- 
speed ratio (or Mach number) than will be obtained in a continuum flow 
field. Since the continuum curves also represent the temperature ratio 
of a gas that has been accelerated adiabatically from rest to any value 
of molecular speed ratio, it may be inferred that, an insulated cylinder 
located in a low— density wind tunnel under free— molecule— flow conditions 
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will attain a temperature which exceeds the total temperature of the gas 
stream. This result, of course, is in direct contrast to the correspond- 
ing phenomenon which occurs under continuum flow conditions where an 
insulated body can have, at most, a temperature equal to the stream stag- 
nation temperature and normally does not even attain this temperature due 
to outward heat flow in the boundary layer. The anomaly can be explained, 
however, by a consideration of the magnitudes of the incident and re— • 
emitted molecular energy. The incident molecular energy is computed 
using the total velocity resulting from combination of the stream mass 
velocity and the random thermal velocities. When the total— velocity 
term is squared there results a cross-product term, 2UV. This cross- 
product term effectively increases the apparent internal energy of the 
gas from the continuum value of (3/2)kT per molecule to a value from 
2kT to (5/2)kT, depending upon the speed and orientation of the body. 

The apparent internal energy becomes (5/2)kT for high speeds and body 
angles of attack greater than zero. This is just equal to the interal 
plus potential energy per molecule of a cube of continuum gas. Therefore, 
the energy incident on .the body in free-molecule flow becomes equal to 
that of a continuum for large values of the molecular— speed ratio for 
surfaces inclined at angles of attack greater than zero. The molecular 
energy which is re— emitted from the surface is assumed to be equal to 
the. energy of a Maxwellian stream issuing effusively from a gas, in equi- 
librium at some yet unspecified temperature, into a perfect vacuum. For 
the case of • an insulated body the temperature of this gas is that of the 
body. The energy of the re-emitted stream calculated in this manner is 
given by equation (A12) to be equal to 2kT w per molecule. The corre- 
sponding energy for a continuum gas at the same temperature is (5/2)kT w 
per molecule if both thermal and potential energies are considered. 

Because of the net energy deficiency of the re-emitted effusive stream 
for a -given wall temperature, a body in free-molecule flow will attain 
a temperature 1 higher than that attained in a continuum flow. 

It is seen from figure 18 that, theoretically, a higher cylinder 
temperature is predicted for the case of a monatomic gas than for a 
diatomic gas. The explanation lies in the additional energy which the 
diatomic gas is able to remove from the cylinder by virtue of the inter- 
nal energy component due to molecular rotation. Of course, the impinging 
diatomic gas transmits more energy to the cylinder due to the rotational 
energy component; however, the rotational component is a smaller fraction 
of the total impinging energy than Of the total re-emitted energy. 

An ex amin ation of the equations for the drag coefficient of an 
insulated cylinder (equations (17) and (l8)) shows that the total drag 
coefficient, Cp = Cp^ + Cp r , is a function only of the moleculai^speed 
ratio. A plot of the total drag coefficient as a function of the molec- 
ular-speed ratio for an insulated cylinder is shown in figure 19- It is 
notable that the drag coefficient approaches infinity as the molecular- 
speed ratio approaches zero. Hence, if the expression for total drag is 
written in the conventional drag-coefficient form, it can be shown from 
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equations (17) and (l8) that Cp becomes inversely proportional to the 
speed ratio as s approaches zero. It is also notable that the drag 
coefficient approaches a constant value as the speed ratio becomes large. 
This means , sinply, that -the effect of the random thermal motions becomes' 
negligible as the speed of the cylinder relative to the gas becomes high. 


Heat-Transfer Test Results 


The salient point of the analysis concerning heat transfer. to a - 
cylinder in a free-molecule field was the prediction that an insulated 
model would attain a temperature in excess of the tunnel stagnation 
temperature. This result was observed in all runs. The maximum rise, 
of the model center-point temperature above the tunnel stagnation 
temperature was 65° F, using nitrogen as the test gas, and 147° F, using 
helium as the test gas. 

' The results of the heat— transfer tests are shown in figures 20 and 
21 in which the ratio of the cylinder— center temperature to the free— 
stream static temperature is plottetl as a function of the molecular— speed 
Iratio for the nitrogen and helium runs, respectively. The stream— static 
temperature was computed from the adiabatic relationship , using the 
measured value of Mach number (or molecular— speed ratio) corresponding to 
the nozzle pressure ratio as determined from the nozzle— calibration tests. 
The temperature ratio curve for continuum flow assuming a temperature- 
recovery factor of unity is shown for comparison. It will be noted that 
in all cases the temperature ratio exceeded the maximum which could be 
attained under continuum flow conditions. It can be seen that, in 
general, the test points indicate that the cylinder did not attain .the 
full temperature rise that was predicted by the free— molecule analysis. 
This result may be ascribed to heat losses occurring from radiation and 
from conduction along the wire to the end supports. The fact that a 
large temperature gradient existed along the length* of the wire due to 
end losses is clearly shown by the cylinder— end. temperatures which are 
listed, together with the cylinder center-point temperature, in table II. 


Drag Test Results 

The measured drag force on the cylinder was reduced to conventional 
dimensionless drag-coefficient form by means of the equation 

Cl) " 7 LrpM 2 ” 2Lrps 2 ( 20 ) 

The data points obtained from equation (20) are shown in figure 22. 
For comparison, the theoretical curves of drag coefficient as a function 
of molecular— speed ratio for an insulated cylinder, taken from figure 19 , 
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are shown as a solid line. It can "be seen that agreement between theory 
and experiment is good considering the limitations of the theory and the 
difficulties attendant with the measurement of drag forces from 1 to 20 
milligrams. The experimental values for the drag coefficient tend to 
fall below the theoretical values at molecular— speed ratios above about 
1.0 and to be higher than the theoretical values at speed ratios below 
1.0. The former characteristic may be due to a tendency toward specular 
reflection at the higher speed ratios while the latter characteristic, at 
subsonic values of speed ratio, may have been caused by interference with 
the flow caused by the shields at the ends of the model. 

The major premise of the existence of free-'inolecule flow was borne 
out by several factors in addition to the good agreement between the 
cal ciliated and experimental values of drag coefficient. One of the 
basic assumptions of the drag analysis is that the drag coefficients 
are independent of Reynolds number or Knudsen number. In contrast, the 
continuum and slip-flow regimes are characterized by the dependence of 
drag coefficients both upon Reynolds number and Mach number (cf., 
reference 16 ). Although it was not experimentally possible to ■ vary the 
Knudsen number independently of the speed ratio for one test gas, the 
use of both nitrogen and helium gave overlapping values of speed ratio 
with an approximate 2:1 variation of Knudsen number. The variation of 
IQiudsen number covered by the tests with molecular-speed ratio is shown 
in figure 23. No trend is noticeable in the drag coefficients obtained 
in each case; consequently, it can be concluded that f re e-molecule flow 
existed. 

A comparison between measured drag coefficient and the drag coeffi- 
cient calculated by methods as given by Lamb (reference 17) illustrates 
the inapplicability of continuum theory in the free-molecule regime. A 
sample calculation made, by Lamb's continuum method for conditions of 
molecular speed ratio equal to 0.5 and a Knudsen number of ' 185 (based on 
cylinder radius) gave a drag coefficient of 390. The corresponding 
experimental value was l4.6. 


CONCLUSIONS 


The following conclusions may be drawn from a comparison of the 
test results with free- molecule-flow theory shown in this report: 

1. The salient point of the heat— transfer analysis, which was the 
prediction that an insulated cylinder would attain a temperature higher 
than the stagnation temperature of the stream, was confirmed by the 
test results. 

2. Drag coefficients, calculated from free-molecule-flow theory, 
agreed well with the measured drag coefficients. As predicted 
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theoretically, the measured drag coefficients were independent of the 
Reynolds and Knudsen numbers and depended primarily upon the molecular- 
speed ratio . 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., Sept. 12, 1950- 
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APPENDIX A 


THE MOMENTUM AND ENERGY' EXCHANGE PROCESSES ACCOMPANYING 
MOLECULAR EMISSION FROM A SURFACE 


For the case of diffuse scattering of a stream of molecules from a 
surface, all directional history of the impinging stream is erased upon 
impact with the surface. The impinging stream is assumed to he scattered 
in a random manner which will he described in detail later. The energy 
and momentum carried hy the scattered stream depends upon the velocity 
of emission which can be related to the energy (or temperature ) level 
that the stream has. The effective temperature or energy of the scat- 
tered stream depends upon the efficiency of the energy— exchange process 
that occurs between the solid surface and the impinging stream. This 
efficiency can be described by introduction of the accommodation coef- 
ficient which is defined as 


a 



• E i “Ew 


(Al) 


where the symbols have the following definitions: 

E^ rate of incident molecular energy 
Ep rate of re— emitted molecular energy 

Ew rate of re-emitted molecular energy that would be carried by the 

scattered stream if it were a stream issuing from a gas in equilib- 
rium at the surface temperature, T w 

It is assumed that the emitted stream has a Maxwellian distribution 
of speed corresponding to a gas in equilibrium at an unspecified tempera^ 
ture, T r . This fictitious gas has, also, an unspecified number of mole- 
cules per unit volume, N r . Now, as shown in reference 8, page 62, the ' 
number of molecules scattered from unit surface area per second that 
have speeds in the range dv and that are moving in a direction that 
makes an angle lying in the range dG with the nor mal to the surface 
is given by 


dn^. = 2itN r A r v 3 e“^ : 


2 2 


sin G cos G dv dG 


(A2) 


Then, the normal component of force imparted to the surface by this group 
of molecules is * ' 


n 2 2 

4 “Px* v 2 

dGr = 2«N r A r mv e sin G cos' G dv dG dA 


(A3) 
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and the translational energy carried by this group is 


5 -0r 2v2 
dE r = nN^mv e 


sin ft cos ft dv dft cLA 


(A4-). 


By integration of these expressions with respect to ft and v, the follow- 
ing values are obtained for the total momentum and energy carried by the 
scattered stream from dA: 


dG r = 2rtN r A r m dA / v e 


o 2 2 

00 4 Py* V 


dv 


I’ 


«/, 


2 2 N~ja 

sin ft cos ft dft = — ■ „ dA 
’ 4p r 2 


i ^ oo es — Bv. V n Jt/2 NuQl 

dE r = TtNjAjjn dA / v e- dv / sin ft cos ft dft = - — dA 

Jo 


2 ^P. r 3 


(A5) 

A 

(A6) 


In order to obtain a value for N r in terms of known variables, the 
number of molecules scattered from the surface is equated to the number 
of molecules ' striking the surface. The number of molecules per second 
scattered from unit area is 


ny = 2nN r A r 


L 


o 2 2 
00 3 V 

v e r dv 


L 


n/2 


sin ft cos ft dft = 


M r 

2 p 1 


(A?) 


and the number of molecules striking unit area per second is shown in 
reference 4 to be given by 


n 


NX 

2 P 


(A8) 


The value of N r is obtained by equating n_ and n from equations 
(AT) and (A8), thus, 

N r = NX (A 9 ) 

By substitution of this value of N^ in the momentum and energy equa- 
tions, the following is obtained; 


dG r 


NmX 

4p3 T 


dA 


- (A10) 


dEL, = 


NmX 


r ~ 2 33 r 2 


dA 


An identical analysis yields a value for of 

NmX dA 


= 


2 


= 2nkT w dJ 


(All) 


(A12 ) 
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For the case of a diatomic gas where the expression for the re-emission 
energy per molecule contains an additional term having the value kT w , 
which accounts for the rotational energy per molecule, the total re- 
emission energy has the value SnkTydA. 


Equation (A10) describes the normal re-emission pressure on the 
front side of a unit inclined area. In order to obtain the net drag 
force (or pressure) in the direction of the mass velocity U, the pres- 
sure on the rear side of the inclined area must be taken into account. 
The computation of the pressure on the rear side of the area is made in 
a similar mann er to the computation of the pressure on the front side of 
the area and leads to an expression 


dG r ' 


NmX* 

4pp r 


dA 


(A13) 


Then, the net drag force due to re— emission pressure in the direction of 
U, the mass velocity, is 


dG r - dG r ‘ = - Q - (X-X 1 ) dA (Al4) 

Replacing X and X* by their values, equation (Al4) becomes 

- *V . ■&. M g^ ln2 6 dA (Al 5 ) 

This expression can be put into drag-coefficient form by use of the 
following substitutions: Nm = p, s/p = U. Equation (A15) then becomes 

dG r -d<V = (A16) 

where s r is the molecular-speed ratio of the re-remitted stream. Then 
the re— emission drag coefficient, based on the total area of one side of 
the plane area, is , 

% T -- C h <w> 

In order to obtain an expression for s r which contains measurable 
or calculable physical quantities it is necessary to write an equation 
expressing an energy balance on the body in question. For the case of an 
insulated cylinder, if radiant energy exchange is neglected, only molecu- 
lar energies are involved and the incident molecular energy can be equated 
to the re— emitted molecular energy. To write an energy balance equation 
for a transverse cylinder, it is assumed that the energy balance is main- 
tained for the cylinder as a whole and not for individual elements of 
the cylinder. This means, essentially, that the cylinder has a uniform 
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temperature, circumferentially. The total incident molecular energy on 
the cylinder is, for a monatomic gas. 




n/2 


X» (a 2 + ) de 


The total molecular energy emitted from the cylinder surface 3s 

2rLp / r 

‘ 2 U 


®r = 




r */ z 
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’de^ 


(ai8) 
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The values of the integrals are shown in appendix B to he 
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When E^ and Ej. are equated from equations (Al8) and (A19) and the 
values of the integrals are' substituted, the following result is obtained: 
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Therefore, 


P r = P 


2(Zi+Z2 ) 


/Z 1 (s 2 +2) 


+ Z 2 ^s 2 + 


(A21) 


A similar analysis for the case of a diatomic gas where the incident 
molecular energy is 


E. = 


r - p - 3 r (s 2 +i) r w ^xd0+(s 2 +i) r R/ x*d0+ r / 
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28 


NACA TN 2244 


and the re— emission energy is 



gives a value of 0 r of 



3(z 1+ z 2 ) 


y 2 i(s 2 +3)+z 2 



(A23) 


(A24) 


Now, the total re-emission drag force acting on dA has been shown 
(cf., equation (Al6)) to he 


dG r - dG r ' . 0 ^ 

To obtain the total re-emission drag on a cylinder, equation (Al6) must 
he integrated over the cylinder surface, thus. 


(arL) vX 


/r r* ,s 


sin 2 6 d0 
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(A25) 


Then, the re-emission drag coefficient based on the projected cylin- 
der area. is 

3/2 

(A26) 


^r ‘ 4Z 


The final equation for the re-emission drag coefficient for an insulated 
cylinder in a monatomic gas stream is obtained by using the relation 
s r = UP r and taking P r from equations (A21) and (A24). 
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and, for an insulated cylinder in a diatomic gas stream, as 



(A2 8) 
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APPENDIX B 


HEAT TRANSFER TO A TRANSVERSE CIRCULAR CYLINDER 
IN FREE-MOLECULE FLOW 


It has been shown (see Analysis) that an energy balance on a cylin- 
der located in a stream of monatomic gas flowing with a velocity U per- 
pendicular to the cylinder axis is expressed by 


2arL 


r /mO 2 \ r rt /s /mU 2 \ 

J n^— + 1 (fkTj d e+y n« A ~ + *«kT \ d0 


If-arLkT- 



Jt/2. p n / s \ 

' nd0 + n’d0 i 


2JfrL [Q-ea(T c 4 -T t 4 )] = 0 


(Bl) 


Equation (Bl) can be simplified by dividing through by 2«rL and replac- 
ing n and n* by their values from equations (2) and (4), thus. 


(mU 2 — 4kT c ) ^ Xd0 + J* 1 * X‘d0^ 

rt/2 


'Jt/2 


'*/ 2 


akTNv, 


2« 


r 2 p Jt/j 

372“ / + 

o ~o 



y' n/fe X»jf» d0^ -eo (T c 4 -T t 4 ) + Q = 0 

n ' 


(B2) 


The integrals to be evaluated are 

r n/ 2 

^o 

jt/2 


jt/a 

u o 

r 

I v 


Xd0 ’ 

(B3) 

X’d0 

(B4) 

Xtyd0 

(B5) 

X*t|f»d0 

(B6) 
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If the value of X in equation (B3) is replaced by its definition (from., 
list of symbols), equation (B3) becomes 



f — a 2 sin 2 


0 


+ 8 8 i n 9 + s sin 6 erf (s sin 9) 


d9 (by) 


The first term of this integral can be evaluated by setting 

s = 2x and cos 29 = z 


Then, d0 


dz 

2 sin 29 


and the limits become 


V 


z = 1 when 0=0 
z = — 1 when 0 = — 


Thus we can write 


n/s 


-s 2 sin 2 0 


e 


d0 


e-x p 1 e xz 

= 2 4 


dz 


(b8) 


Wow, by definition (cf., equation (11.), p. 46, reference 18), the modified 
Bessel function of order n can be written as 


[ n (x) = — — - ( - YT 1 e±xz (l-z 2 ) 

r[n+(l/2)] V 2 yj 


n-(i/2) 


dz 


Then 


i 0 u) = — - — r 1 


y»r(i/ 2 ) J 


dz 


and since 


I) 


(B9) 


( 
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the first term, of integral (B3) becomes 


n P -s 2 /? T ( s 2 

2 ® 10 U 


' (BIO) 


The second term (cf . , equation (B7)) is 


s / ' sin 0 d0 = a 

J o 


jr/s 


The third "term of equation (B7) is 


Jt/2 


*/it~ s / ' sin 0 erf (s sin 0) d0 

J o 


which can be integrated by parts. Let 


and 


Then 


and 


erf (s sin 6 ) = u 


sin 9 d0 = dv 


2s -e 2 sin 2 0 
du = - — e cos 0 d0 

nr 


v = -cos 0 


Then, % the integral becomes. 


(Bll) 


n/n 


— cob 9 erf (s sin 0) 


-.it/a 


+ 2s‘ 

-*o °o 


rt/2 


cos 2 0 e - ® 2 sin20 d0 


(B12) 

t 

The bracketed term of equation (B12) is zero, leaving the second term 
yet to be evaluated.. Make the substitutions 
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cos 2 a = it££ 2 Ji 
2 

sin 2 0 = i=52L Ji 
2 

s 2 = 2x 


Then 


48 -f 


and., replacing limits 


0=0 


0 =| 


u=0 

u=« 


there results 


2 P rt / 2 2 ~ -e 2 sin 2 0 s 2 P 11 , 

a J cos 2 0 e d.0 = J (l+cos u) e 


— x( l-cos u) 


du 


s 2 e r — x cos u, s 2 e X / x cos u , 

e du+ — ^ — / cos u e du 


2 J Q - — 2 


I 


(B13) 


Integral (B13) becomes, when substituting z=cos u and changing limits, 

s 2 e~ x P 1 e xz 
2 


-1 / 1 hs 2 


a 2 -x , /•■" XZ 

dz + LL_ / £2 dz 


«sf. -fl 2 ^ 
2 e 


2 - 

( 4 ) + ^ ( 4 ) 


(Blh) 


which is obtained -by integrating the second term of (Blh) -by carts and 
substituting equation (B9). Then, the final value for integral (B3) is 


n n / 2 « -s 2 /2 

j Xd0 = -e 


T ( s2 '^ /T «S 2 — S 2 / 2 _ f B 2 \ ( g 2 ^ 

Z o {- 2 ) + > 8 + — ’ e L J o W + 11 V2V 


(B15) 
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Now, integrals (b 4), (B5) > and (B6) contain only terms which can be eval- 
uated by the methods shown above by which integral (B3) was evaluated. 

The results are 


r 


X*de = 


Jt -«2 
2 6 


■/• (f ) [i 0 ($) ♦ ii (f); 


J xtae-« tT**! 2 i 0 (^) + |-/7 s+ ^ « s^ - ® 2 / 2 i G (^^) 


■jC 


n/2 


X»i|r*de=Jte' 




72 ^(f) - \-P s+ f - 2 ^ 2/2 [lo(# * 


(B16) 


+ Ii K-o 


(BIT) 


(B18) 


In order to eliminate an excessive amount of writing, new variables are 
defined as follows: 


Zi = * e - ® 2 / 2 I 0 (B19) 


Z 2 = at s^ - ® 2 / 2 £ I 0 + Ii ] (B20) 

Values of the dimens ionless variables Z x and Z 2 were confuted 
from thq tables of Bessel functions in reference 19 . Then, equation (B21) 
is obtained by substitution of equations (B15) to (B20) into equation (B2) 

Nv^jci* akTNv^ / c \ 

- ^375 ( inti 2 — 4kT c )(Zi+Z 2 ) + -gyg | Z S J - ea(T c 4 -T t 4 ) + Q=0 (B21) 


If mv m 2 /2 = KT, KTN = p, and s = U/v m are substituted into equation 
(B2l), the following equation results: 
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2 ^ (Z x +Z 2 ) 




Zi(s 2 +2) + Z2 


( s2+ ?)] 


+ WB- [«<’(T c 4 -v>- <i ) =o < bs2 > 


which' is the final energy exchange equation for a monatomic gas. An iden- 
tical treatment of equation (13) for a diatomic gas gives 

3 ^ (Zi-tZs)- [z x (s2+3)+z 2 ( s 2 + 0J + [ 6a (T c ^ t 4 ) - Q] = 0 

The groups f(s)=Z x (s 2 +3) + Z 2 (^s 2 + 0 , g(s)=3(Z x +Z 2 ), f°(s)=Z x (s 2 +2) + 

Z 2 (" s2+ D , and g°(s)=2(Z x +Z 2 ) are tabulated in the following table: 


s 

f(s) 

g(s) 

f°(s) 

g°(s) 

0.2 

9.80271 

9.61234 

6.59860 

6.40823 

.4 

10.94784 

10.16410 

7.55981 

6.77607 

.6 

12.89399 

11.04924 

9.21092 

7.36616 

.8 

15.69746 

12.22282 

11.62319 

8.14855 

1.0 

19.43609 

13.63288 

14.89180 

9.08859 

1.2 

24.20781 

. 15.22758 

19.13179 

10.15172 

1.4 

30.12863 

16.96068 

24.47507 

11.30712 

1.6 , 

37.32923 

18.79408 

31.06454 

12.52939 

1.8 

45.95230 

20.69872 

39.05273 

13.79915 

2.0 

56.14932 

22.65358 

48.59807 

15.10239 

2.2 

68.07771 

24.64392 

59.86307 

16.42928 

2.4 

81.89947 

26.65990 

73 .01284 

17.77327 

2.6 

97.77902 

28.69470 

88.21412 

19.12980 

2.8 

115.88287 

30.74382 

105.63493 

20.49588 

3.0 

136.37821 

32.80398 

125.44355 

21.86932 

3.2 

159.43361 

34.87300 

147.80928 

23.24867 

3.4 

185.21769 

36.94914 

172.90131 

24.63276 

3.6 

213.89927 

39.03112 

200.88890 

26.02075 

3.8 

245.64801 

41.11800 

231.94201 

27.41200 

4.0 

280.63340 

43.20904 

266.23039 

28.80603 

4.2 

319.02474 

45.30354 

303.92365 

30.20236 

4.4 

360.99204 

47.40108 

345.19168 

31.60072 

4.6 

406.70528 

49.50126 

390.20486 

33.00084 

4.8 

456.33364 

51.60366 

439.13242 

34.40244 

5.0 

510.04756 

53.70808 

492.14487 

35.80539 

5.2 

568.01730 

55.81428 

549.41254 

37.20952 

5.4 

630.41182 

57.92200 

611.10449 

38.61467 

5.6 

697.40229 

60.03114 

677.39191 

40.02076 
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Then, integration of equation (Cl) and. substitution of the relations of 
equation (C2) yield 


7LG± = £*f. 


2 e 

L Jt 1 / 2 


— s 2 sin 2 0 


+ 2 sin 0 erf (s sin 0) 


(i *&1 


dA (C3) 


It is remarkable that a purely analytic solution for this type flow leads 
to a simple equation of the type of (C3) in which the bracketed term is a 
dimensionless function of molecular speed ratio and which 1 may be inter- 
preted as an ordinary drag coeffient. 

In order to obtain the drag force on a transverse cylinder due to 
impinging momentum, equation (C3) ,- which was derived for an elemental 
'plane area, must be integrated over the surface of the cylinder. To this 
end, dA is written as 


dA = rid© 

Then, for the cylinder, the impinging drag force is 


(04) 


G^ = 2rL 




J/ s r 2e -e z Sin 2 0. 


- 1/2 
it s 


+ 2 sin 0 erf(s sin 0) (l + — i— j 

v 2s 2 / 


d0 
(C5) 

This integral can be evaluated by the methods shown in appendix B. The 
final value for the drag force due to impinging molecules is 


= 2rL 


[■•(?)*'.(?)]} >*> 


Thus, there results from equation (C6) an impinging drag coefficient 
based on the projected area of the cylinder, which has the value 


C Di = 


l/2 

jt < e 


-s 2 /2 




(C7) 
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TABLE II.- TEST DATA 



S ® tSf 

■S ^ m 

CQ CQ ri 
_ W O 
§ <D ^ 

3 u o 

<D ft -H 

S 1 

ca 


Pl ^ 
o q) . bo 

33 

ti CQ <H 

is ° 

.a} U g 
-P ft § 



N 2 

2.49 

71.0 

72.8 

2.70 

73.0 

72.6 

101.3 

127.7 

n 2 

2.78 

120.0 

120.2 

5.00 

73.0 

72.2 

94.0 

119.0 

He 

.55 

18.0 

19.5 

.70, 

70.2 

70.0 

80.0 

87.0 

He 

2.08 

51.0 

49.4 

2.15 

72.2 

70.2 

123.7 

189.0 

He 

2.37 

89.0 

95.6 

4.00 

64.0 

62.2 

120.0 

211.0 

n 2 

2.07 

71.6 

74.4 

1.20 

75.5 

75.2 

103.3 

131.7 

n 2 

2.25 

119.0 

132.7 

2.20 

74.7 

75.2 

95.7 

122.7 

He 

1.57 

53.5 

52.0 

1.00 

75.5 

75.5 

105.7 

141.5 

He 

1.75 

91.0 

91.0 

1.70 

65.0 

67.5 

114.4 . 

183.3 

n 2 

1.85 

68.7 

68.7 

.90 

75.0 

74.4 

101.7 

129.5 

He 

.73 

32.0 

33.0 

.35 

76.2 

74.0 

82.5 

93.7 

He 

1.00 

44.3 

46.5 

.60 

67.4 

68.7 

87.6 

110.5 

He 

.89 

22.6 

22.6 

.19 

68.5 

74.7 

80.7 

87.5 

He 

1.67 

55.0 

55.0 

.50 

74.0 

74.4 

101.0 

132.0 

He 

1.68 

78.5 

78.5 

■ .70 

68,2 

69.4 

107.0 

162.5 

He 

1.73 

102.0 

102.1 

.90 

70.7 

70.8 

110.7 

174.5 

He 

1.79 

117.0 

122.1 

1.05 

71.2 

71.5 

108.7 

175.7 

He 

1.83 

130.0 

135.0 

1.15 

70.8 

71.2 

106.3 

173.5 

H? 

.78 

18.6 

18.6 

.15 

73.7 

74.3 

76.6 

80.5 

He 

1.23 

43.5 

46.0 

.30 

76.3 

74.3 

94.2 

117.7 

He 

1.36 

63.6 

62.5 

.40 

74.2 

74.3 

102.7 

140.3 

He 

1.44 

80.5 

80.5 

.47 

73.3 

74.8 

105.0 

147.8 


Average Mach number over center half inch of stream. 
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Figure 3.- impact- pressure- measuring system. 
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Figure 4- Static- pressure-measuring system. 



NACA TN-2244 


Note: All dimensions given in inches 
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Time, sec 


Figure 6 -Time response of the static- pressure -measuring system. 
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Figure 7 .— Micro-photograph of model showing weld at junction of iron and 
constantan. Enlarged 60 times. 








Adjustable weights 


NACA TN 22kk 


53 



Figure 9- Diagrammatic sketch of drag balance. 
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Distance from center of stream, inches 


Figure //. - Variation of impact and static pressure with 
nozzle radius for nozzle number / using nitrogen gas. 
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Figure 13.- Comparison of Mach number distribution in nozzle 
number / as calculated with and without viscous correction, 
test gas, nitrogen . 
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Figure 14.- Comparison of Mach number distribution in nozzle 
number / as calculated with and without viscous correction , 
test gas, helium . 
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Figure 15 *— Shock wave from 5 ° half-angle wedge at pressure of Ilk microns 
as seen with nitrogen glow discharge. 
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0 .4 .8 12 1.6 2.0 2.4 

Molecular speed ratio , s 

Figure 18- Variation of the ratio of cylinder temperature to free - 
stream static temperature, with molecular speed ratio for an 
adiabatic cylinder in free -molecule flow. 
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Molecular speed ratio, s 

Figure 19 - The total drag coefficient for an insulated transverse 
cylinder in a free- molecule flow as a function of the molecular 
speed ratio, s. 
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Molecular speed ratio, s 

Figure 20 - Variation of the ratio of model temperature to 
free - stream static temperature with molecular speed ratio 
for a diatomic gas. 
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0 .4 .6 1.2 !6 2.0 

Molecular speed ratio, s 

Figure 21- Variation of the ratio of mode! temperature to free - 
stream static temperature with molecular speed ratio for a 
monatomic gas . 
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Molecular speed ratio, $ 


Figure 22- Variation of total drag coefficient with molecular 
speed ratio. 
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Figure 23- The range of Knudsen numbers and molecu- 
lar speed ratios covered in the tests. Knudsen 
number based upon model radius. 
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